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Using gauge/gravity duality, we study the creation and evolution of boost invariant anisotropic, 
strongly coupled J\f = 4 supersymmetric Yang-Mills plasma. In the dual gravitational description, 
this corresponds to horizon formation in a geometry driven to be anisotropic by a time-dependent 
change in boundary conditions. 



I. INTRODUCTION 

The study of non-equilibrium phenomena in QCD and 
other non-Abelian quantum field theories is a topic of 
much interest, with applications to heavy-ion collisions, 
early universe cosmology, and other areas. Much has 
been learned about near-equilibrium dynamics at weak 
coupling where a quasiparticle picture is valid Q3 [2] 
IU IH |5], and there has also been considerable progress 
understanding certain theories at very strong coupling 
[HI B HI El] j thanks to the development of gauge/gravity 
duality jTUJ, [TTJ Q2] . But very little progress has been 
made in regimes where a theory is both strongly coupled 
and far from equilibrium. 

Heavy ion collisions at the Relativistic Heavy Ion 
Collider (RHIC) are believed to produce a deconfmed, 
strongly coupled quark-gluon plasma (QGP) [131 03] ■ In 
the initial stages of the collision, during which the QGP 
is produced, the system is surely far from equilibrium and 
cannot be described by hydrodynamics. However, mod- 
eling based on near-ideal hydrodynamics strongly sug- 
gests that a hydrodynamic treatment becomes applica- 
ble rather quickly, perhaps on times < 1 fm/c after the 
collision event |15j . Understanding the dynamics respon- 
sible for such a rapid approach to local equilibrium from 
a far-from-equilibrium initial state is a challenge. 

Colliding nuclei at sufficiently high energy is the only 
experimentally accessible approach for creating quark- 
gluon plasma. However, analyzing the dynamics — from 
the creation of the initial highly non-equilibrium state, 
through partial equilibration, hydrodynamic evolution, 
hadronization, and eventual freeze-out — via a first prin- 
ciples calculation in QCD is not currently possible. Nev- 
ertheless, aspects of this process involving strongly cou- 
pled dynamics can be studied in a controlled setting in 
a class of theories which describe non-Abelian plasmas 
similar to the QGP, and which possess dual gravitational 
descriptions. The best known example is M = 4 su- 
persymmetric Yang-Mills (SYM) theory [TO] . In this 
theory one can study the collision of shock waves using 
gauge/gravity duality (THJ, [XT] [T5J QI5]. The shock waves 
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have a very small thickness along the collision axis and 
can be localized in the transverse directions [20] . There- 
fore, qualitatively at least, they resemble the Lorentz- 
contracted relativistic nuclei in a heavy ion collision. 

In the dual gravitational description, collisions of shock 
waves in SYM turn into a problem of colliding gravita- 
tional shock waves in five dimensions. The resulting bd 
numerical relativity problem is still quite challenging, but 
may be feasible using techniques which are adapted from 
current work in 4c? numerical relativity. One purpose of 
this paper is to begin exploring some of the needed adap- 
tations, albeit in a setting which is simpler than colliding 
shock waves. 

The immediate goal of this paper is to study how 
quickly a far-from-equilibrium strongly-coupled non- 
Abelian plasma relaxes to a regime in which a hydrody- 
namic description is accurate. The answer to this ques- 
tion will necessarily have some sensitivity to how the 
initial state is created. A conceptually simple way to 
prepare non-equilibrium states is to start in the ground 
state, and then to turn on time-dependent background 
fields coupled to operators of interest. After the back- 
ground fields are turned off, one can then watch the 
subsequent evolution of the system. Since a hydrody- 
namic description requires that the stress tensor (in the 
local fluid rest frame) be nearly isotropic [3], particu- 
larly interesting initial states are those in which the ini- 
tial stress tensor is driven to be highly anisotropic. A 
natural way to do this is to make the spatial geometry 
in which the field theory lives be time dependent and 
anisotropic [2"T1 |2"2"] . 

At weak coupling, the addition of energy to the ground 
state by a time-dependent gravitational field can be 
understood in terms of particle production. A time- 
dependent spacetime geometry will create quanta [23], 
and if the time dependence of the deformation in the ge- 
ometry is anisotropic, then the momentum distribution 
of created quanta will be anisotropic as well. After the 
geometry ceases to evolve, quanta will continue to col- 
lide and interact and eventually (on a timescale which at 
weak coupling diverges like 1/A 2 , with A the 't Hooft cou- 
pling) the system may approach a state in approximate 
local thermal equilibrium. 

This quasi-particle picture breaks down as the strength 
of the coupling increases, and one must understand 
the process of plasma production and relaxation using 
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a different physical description. For large N c SYM, 
gauge/gravity duality provides an alternative picture in- 
volving black hole formation in five dimensions. As we 
discuss in Section [TTJ the gravitational dual will involve 
a 5d curved spacetime with a 4c? boundary which has a 
time dependent geometry. The boundary geometry cor- 
responds to the spacetime geometry of the SYM field 
theory. A time-dependent deformation in the Ad bound- 
ary geometry will produce gravitational radiation which 
propagates into the fifth dimension. This radiation will 
necessarily produce a black hole [5T]. It is natural that 
the gravitational description of plasma formation and re- 
laxation involves horizon formation, since at late times 
the system will be in a near-equilibrium state with non- 
zero entropy. 

The presence of a black hole acts as an absorber of 
gravitational radiation and therefore, after the produc- 
tion of gravitational radiation on the boundary ceases, 
the 5d geometry will relax onto a smooth and slowly 
varying form. This relaxation is dual to the relaxation 
of non-hydrodynamic degrees of freedom in the quantum 
field theory [9]. Therefore, by studying the evolution of 
the 5d black hole geometry, one can gain insight into the 
creation and relaxation SYM plasma. 

For simplicity, in this paper we limit attention to 4d ge- 
ometries which have two dimensional spatial homogene- 
ity and 0(2) rotation invariance in the x± = {x 1 , x 2 } 
directions, and which are invariant under boosts in the 
xii = x 3 direction. As we discuss in Section Nlj this 
reduces the gravitational dynamics to a system of two- 
dimensional PDEs, which we solve numerically. Besides 
making the gravitational calculation simpler, these as- 
sumptions serve an additional purpose. With these sym- 
metries, the late time asymptotics of the 5c? geometry 
(and the corresponding asymptotics of the stress tensor) 
are known analytically [231 HS1 US] • We will therefore be 
able to compare directly our numerical results, valid at 
all times, to the known late time asymptotics. 

Boost invariance implies that the natural coordinates 
to use are proper time r and rapidity y (with x° = 
t cosh?/ and x\\ = t sinhy). In these coordinates, the 
metric of Ad Minkowski space (in the interior of the r = 
cone) is ds 2 = —dT 2 + dx\ + T 2 dy 2 . A deformation of the 
geometry, respecting the above symmetry constraints, in- 
duced by a time-dependent shear may be written in the 
form 

ds 2 = _ dT 2 + e7 (r) dx 2 ± + T 2 e -2 7 (r) dy 2 ^ 

The function 7(7") characterizes the time-dependent 
shear; neglecting 4d gravity, 7(7-) is a function one is 
free to choose arbitrarily. For this study, we chose 

7 (r) = cG (1 - (t-t ) 2 /A 2 ) [1 - (r-r ) 2 /A 2 ] 6 

xe -i/[i-(T-r„)VA 2 ] ; (2) 

with the unit step function. (Inclusion of the [1 — 
(t— tq) 2 /A 2 ] 6 factor makes the first few derivatives of 




FIG. 1: A spacetime diagram depicting several stages of the 
evolution of the field theory state in response to the changing 
spatial geometry. At proper time r = Ti, the 4d spacetime ge- 
ometry starts to deform. The region of spacetime where the 
geometry undergoes time-dependent deformation is shown as 
the red region, labeled I. After proper time t = 77, the de- 
formation in 4d spacetime geometry turns off and the field 
theory state is out of equilibrium. From proper time 77 to 
T», shown as the yellow region labeled II, the system is sig- 
nificantly anisotropic and not yet close to local equilibrium. 
After time r*, shown in green and labeled III, the system is 
close to local equilibrium and the evolution of the stress tensor 
is well-described by hydrodynamics. 



7(7") better behaved as t— tq — > ±A.) The function 7(7-) 
has compact support and is infinitely differentiable; 7(7") 
and all its derivatives vanish at the endpoints of the inter- 
val (Tj, 77) , with T{ = t — A and 77 = t + A. We choose 
To = § A so the geometry is flat at r = 0. 1 We choose to 
measure all dimcnsionful quantities in units where A = 1 
(so n = 1/4 and 77 = 9/4). 

Fig. [T] shows a spacetime diagram schematically de- 
picting several stages in the evolution of the SYM state. 
Hyperbola inside the forward lightcone are constant r 
surfaces. Prior to r = 77 the system is in the ground 
state. The region of spacetime where the geometry is 
deformed from flat space is shown as the red region la- 
beled I in Fig. [l] At coordinate time t = Ti the geometry 
of spacetime begins to deform in the vicinity of x\\ = 0. 
As time progresses, the deformation splits into two local- 
ized regions centered about Xn ~ ±t, which subsequently 
separate and move in the ±xy directions at the speeds 
asymptotically approaching the speed of light. After the 
"pulse" of spacetime deformation passes, the system will 
be left in an excited, anisotropic, non-equilibrium state. 
That is, the deformation in the geometry will have done 
work on the field theory state. This region, labeled II, 



Choosing to > A is convenient for numerics as our coordinate 
system becomes singular on the r = lightcone. The particular 
choice to = I A was made so that our numerical results (which 
begin at r = 0) contain a small interval of unmodified geometry 
before the deformation turns on. For an interesting discussion of 
non-equilibrium boost invariant states near r = see Ref. |27| . 
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is shown in yellow in Fig. [T] It is in this region that we 
can study the relaxation of a far-from-equilibrium non- 
equilibrium state. After some amount of proper time t*, 
the system will have relaxed to a point where a hydrody- 
namic description of the continuing evolution is accurate. 
This final hydrodynamic regime is shown schematically 
in green, and labeled III, in Fig. [T] As the late time 
hydrodynamic solution to boost invariant flow is known 
analytically, we choose to define r* precisely as the time 
beyond which the stress tensor coincides with the hydro- 
dynamic approximation to within 10%. 

Our task then is to find and in particular, see how 
it correlates with quantities such as the effective temper- 
ature T* at time r*. In the c —> oo limit, which corre- 
sponds to a diverging size of the deformation in the Ad 
geometry, it is inevitable that approaches r/. This 
is natural in conformal theories, since relaxation times 
of non-hydrodynamic degrees of freedom are set by the 
local energy density, and this diverges when c — > oo. In 
other words, in the limit where 1/T* <C A, the system re- 
sponds adiabatically to the deformation in the geometry 
and non-hydrodynamic degrees of freedom can remain in 
equilibrium. A hydrodynamic description (without driv- 
ing terms) will be accurate the moment the geometry 
stops changing. Hence, in this limit one learns noth- 
ing about the dynamics associated with the relaxation of 
non-hydrodynamic modes. 

More interesting is the case where the effective temper- 
ature satisfies 1 /T* > A. This is the regime we will study. 
Within this regime, the system can be significantly out- 
of-equilibrium after the 4c? geometry becomes flat. When 
this is the case, we find that the entire process of plasma 
creation and relaxation to approximate local equilibrium 
(i.e., to a point where subsequent evolution is accurately 
described by viscous hydrodynamics) occurs over a time 
which varies between one and two times 1/T*. 

This result is consistent with the findings in our ear- 
lier work |21j where we studied isotropization in a homo- 
geneous strongly coupled AA = 4 SYM plasma. In that 
work, all spatial gradients vanished. There was no exci- 
tation whatsoever of hydrodynamic degrees of freedom, 
and the system relaxed exponentially toward equilibrium. 
In contrast, the dynamics of the boost-invariant plasma 
in the present work involves both hydrodynamic and 
non-hydrodynamic degrees of freedom. The results we 
present display a rather clear separation between far- 
from-equilibrium response, which cannot be described 
by hydrodynamics, followed by later "near local equi- 
librium" dynamics which is accurately described by vis- 
cous hydrodynamics. A noteworthy finding is that the 
domain of utility of hydrodynamics is not limited by 
when higher order terms in the hydrodynamic expansion 
become comparable to the lowest order viscous terms, 
rather it is determined by the relative importance of non- 
hydrodynamic degrees of freedom. 



II. GRAVITATIONAL DESCRIPTION 

Gauge/gravity duality 1 0J provides a gravitational de- 
scription of large N c SYM in which the 5c? dual geometry 
is governed by Einstein's equations with a cosmological 
constant. Einstein's equations imply that the boundary 
metric g^ v (%), which characterizes the geometry of the 
spacetime boundary, is dynamically unconstrained. The 
specification of g^ v (x) serves as a boundary condition for 
the 5c? Einstein equations. This reflects the fact that 4e? 
gravitational dynamics is neglected; the dual field the- 
ory residing on the boundary responds to the boundary 
geometry but does not back-react on the 4c? boundary 
geometry. 

Diffeomorphism and spatial 3c? translation invariance, 
together with our assumed 0(2) rotation invariance, al- 
lows one to chose a 5c? bulk metric of the form 



ds 2 =- Adr 2 



X 2 \e B 



dx\ 



,-2-B 



c?y 2 ] + 2dr dr , (3) 



where A, B, and E are all functions of the bulk radial co- 
ordinate r and time r only. The coordinates r and r are 
generalized infalling Eddington-Finkclstcin coordinates. 
Mailing radial null geodesies have constant values of t 
(as well as x± and y). Outgoing radial null geodesies 
satisfy dr/dr = \A. The geometry in the bulk at r > 
corresponds to the causal future of t = on the bound- 
ary. The form of the metric (|3| is invariant under the 
residual diffeomorphism r — > r + f(r), where /(r) is an 
arbitrary function. 

With a metric of the form ([3]), Einstein's equations 
may be written very compactly as 



E(E)' + 2E'E-2E 2 , 
E (B)' + l(E'B + B' £) , 
A" + 3B'B- 12£'E/E 2 
E + | (B 2 E - A' E) , 
E" 



4, 



B' 2 T,. 



where, for any function h(r,r), 



h! = d r h, h = d T h+\Ad r h. 



(4a) 
(4b) 
(4c) 
(4d) 
(4e) 

(5) 



The derivative h! is a directional derivative of h along 
infalling radial null geodesies, while the derivative h is 
the directional derivative of h along outgoing null radial 
geodesies. Eqs. ( 4d I and (4e| are constraint equations; 



the radial derivative of Eq. ( 4d I and the time derivative 



of Eq. (4e) are implied by Eqs. (4a)-(4c 



The above set of differential equations must be solved 
subject to boundary conditions imposed at r = oo. The 
requisite condition is simply that the boundary metric 
g^ v (x) coincide with our choice (Jl| of the 4c? geometry. 
In particular, we must have 



lim £(r,r)/r = T 1/3 , 

r — >oo 

lim B(r, r) = — | lnr + j(r) 



(6a) 
(6b) 
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One may fix the residual diffeomorphism invariance by 
also demanding that 



lim [A(r, r) 



jr = 0. 



(7) 



These boundary conditions, plus initial data satisfying 
the constraint ( 4e I on some r = Tt slice, uniquely specify 



the subsequent evolution of the geometry. 

Near the boundary one may solve Einstein's equations 
with a power series expansion in r. Specifically, A, B and 
E have asymptotic expansions of the form 



A(r,r) = I On (r) + a»(r) log r] r 2 ~", (8a) 

n=0 

B(r,r) = J2[bn(r)+(3 n (T)logr]r- n , (8b) 
S(r,r) = [s«(r)+<7 n (r) log r] r 1 "™ . (8c) 

The boundary conditions d6| and (ITf imply that &o( r ) = 
-§1ht + 7 (t), s (r) = r 1 / 3 , a (r) = 1, oi(t) = 0, 
and that the coefficients of the corresponding logarith- 
mic terms vanish. Substituting the above expansions into 
Einstein's equations and solving the resulting equations 
order by order in r, one finds that there is one undeter- 
mined coefficient, 64 (r). All other coefficients are deter- 
mined by the boundary geometry, Einstein's equations, 
and 64(7-). 2 

Given a solution to Einstein's equations, the SYM 
stress tensor is determined by the near-boundary behav- 
ior of the 5d metric [28 . If Sq denotes the gravitational 
action, then the SYM stress tensor is given by 



SSr. 



(9) 



By substituting the above series expansions into the vari- 
ation of the on-shell gravitational action, one may com- 
pute the expectation value of the stress tensor in terms 
of the expansion coefficients. This procedure has been 
carried out in Ref. |28j, so we simply quote the results. 
In terms of the expansion coefficients, the SYM stress 
tensor reads 



2tt 2 



diag(-£,Vj_,V± t V\\) 



where 



£ = 



1 04 

Z«4 



b 4 + V ± , 
26 4 + % 



(10) 

(11a) 
(lib) 
(11c) 



2 The coefficient 0,4 is determined by a first order ordinary differ- 
ential equation, which can be obtained from the condition that 
the SYM stress tensor be covariantly conserved. All other coeffi- 
cients are determined algebraically from bo(r), &4(t), a^lr) and 
their derivatives. 



and 

£ : 



„->.s ' ■ iTb( 21 7i 2 + 472)T 2 



5 „, „--3 

; 71 r 



jk( 3 7i -87172 -73)7" 1 



96 

5 i f .(37f + 14 72 2 -4 7l73 ), 

fi7ir- 3 -34^(2397 7l 2 
S ig(57 7l 3 + 4887 1 7 2 + 70 7 3)T- 1 
^(21 7 f - 468 7l 2 72 + IO72 + 47173 + 64 74 ) 



(12a) 



1 T - 4 

6 ' 



1444 72 ) t 



'3 ' 



(12b) 

34^(53797? + 2828 72 )r- 2 
2 -|s (120 7l + 458 7l72 + 73 73 ) t- 1 
+ ^(21 7 f + 936 7l 2 72 + 10 7 | + 4 7l73 - 12874) , 

(12c) 



-49 y T -3 
288 '! 

„3 



with 7 „ = d n ^jdi 



III. NUMERICS 



One may solve the Einstein equations (4a|-(4c| for the 
time derivatives E, B, and A". Define 

/>oo 

6(r,r)= / dw[Z(w,T) 3 - hiiw^)] -H^t), 

J r 

(13a) 

(13b) 
(14) 

(15a) 
(15b) 
(15c) 



$(r,r)= / dw 2e(w,T)B'(w,T)Y l (w,T)- 3/2 

J r 

~h 2 (w,T) -H 2 {r,T), 

where H n is an indefinite radial integral of h n , 
h„ = H' . 



Then Eqs. (4a|— (4c| are solved by 



E 
B 
A" 



-26 E , 
-§<I>E~ 3/2 , 
-4 - 240 E'E 



' 4 +|$B'E- 3 / 2 . 



The functions h n (r, r) are not constrained by Einstein's 
equations — their presence inside the integrands of 
Eq. (13 1 are compensated by the subtraction of their 



integrals H n (r,r). However, since we have chosen the 
upper limit of integration in Eq. (13 1 to be r = 00, the 



functions h n (r,T) must be suitably chosen so that the 
integrals (13 1 are convergent. The simplest choice to ac- 



complish this is to set /ii(r, r) equal to the asymptotic 
expansion of E(r, r) 3 up to order 1/r , for some k > 1, 
and to set /i 2 (r, r) equal to the asymptotic expansion of 
2©(r, r)B'{r, r)/E(r, r) 3 / 2 up to order l/r k . In our nu- 
merical solutions reported below, we use k > 4. This 
choice makes the large r contribution to the integrals in 
Eq. (13 1 quite small and consequently reduces cutoff de- 



pendence. As the coefficients of the series expansions fl8| 
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only depend on bo (r) and 64 (r) and their r derivatives, 
this choice determines /in(?"i t) in terms of one unknown 
function 64 (r). 

With the subtraction functions /i n specified by the 
aforementioned asymptotic expansions, integrating Eq. 
( 14 1 fixes the compensating integrals H n up to an in- 



tegration constant which is an arbitrary function of r. 
Integrating Eq. (15c I for A(r, r) introduces two further 
(t dependent) constants of integration. The most direct 
route for fixing these constants of integration is to match 
the large r behavior of the expressions (15a I and (15b) 



and the integrated version of Eq. ( 15c ) to the correspond- 
ing expressions obtained from the series expansions (pj. 
This fixes all integration constants in terms of 60 and 64. 

Our algorithm for solving the initial value problem 
with time dependent boundary conditions is as follows. 
At time Ti the geometry is AdSs with the metric 



ds 1 



-dT 2 + dxl+(T+l) 2 dy 2 



Therefore, at the initial time we have 

B(r,n) = -| M (ri + i) 
E(r,Ti) - r( n +l) 1/3 , 
A(r,Ti) = r 2 . 



2drdr. (16) 

(17a) 
(17b) 
(17c) 



With A(r,Ti), B(r,Ti) and £(r, Tj) known, one can then 
compute the time derivatives d T B(r,Ti) and 9 T E(r, Tj) 
from Eqs. (15b I and ( 15a| , and step forward in time, 



B(r, Ti + At) 
S(r,T, + Ar) 



B(r,Ti) + d T B(r,Ti)AT, (18) 
E(r,r 4 )+a T E(r,r 2 )AT. (19) 



With B(r, r^+Ar) and E(r, Tj+At) known, one can then 
integrate Eq. (15c I to determine A(r, Tj+At). With the 
complete geometry on the time slice r = r^+Ar deter- 
mined, one may then repeat the entire process and take 
another step forward in time. 3 

An important practical matter is fixing the computa- 
tion domain in r — how far into the bulk does one want 
to compute the geometry? If a horizon forms, then one 
may excise the geometry inside the horizon as this re- 
gion is causally disconnected from the geometry outside 
the horizon. Furthermore, one must excise the geometry 
to avoid singularities behind horizons |29j . To perform 
the excision, one first identifies the location of an appar- 
ent horizon (an outermost marginally trapped surface) 
which, if it exists, must lie inside an event horizon 301 . 



3 Because we are working with a discretized version of Einstein's 
equations, the discretized version of the constraint equation 
is not automatically implied by the discretized version of the 
other Einstein equations. To minimize the amount of accumu- 
lated error, we also monitor the accuracy of the constraint equa- 
tion |4a}, and make tiny adjustments to £ to prevent growing 
violation of the constraint. 
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FIG. 2: The congruence of outgoing radial null geodesies. 
The surface coloring displays A/r 2 . Before time Ti = 1/4 
this quantity equals one. The excised region lies inside the 
apparent horizon, which is shown by the dashed magenta line. 
The geodesic shown as a solid blue line is the event horizon; it 
separates geodesies which escape to the boundary from those 
which cannot escape. 



We have chosen to make the incision slightly inside the 
location of the apparent horizon. For the metric |3j, the 
location r^(r) of the apparent horizon is given by th e out - 
ermost point where E(rh(r),r) = or, from Eq. (15a I, 
6Mt),t 



0. 



IV. RESULTS AND DISCUSSION 

We first discuss our results from the 5d gravitational 
perspective and present data for c — 1. Results for 
other values of c are presented below, but the qualita- 
tive features of the results are independent of the value 
of c. Fig. [2] shows a congruence of outgoing radial null 
geodesies for c = 1. The geodesies are obtained by in- 
tegrating dr/dT — ^A(t,t). The colored surface in the 
plot displays the value of A/r 2 . Excised from the plot 
is a region of the geometry behind the apparent horizon, 
whose location is shown by the magenta dotted line. 

At times r < r, = 1/4, the boundary geometry is 
static and A/r 2 = 1. The outgoing geodesic congruence 
at early times therefore satisfies 



+ 2/r = const . 



(20) 



and hence appears as parallel straight lines on the left 
side of Fig. [2] These are just radial geodesies in AdSs, 
which is the geometry dual to the initial zero temperature 
ground state. After time n the boundary geometry starts 
to change, A/r 2 deviates from unity, and the congruence 
departs from the zero temperature form ( |20[ ). 

Perhaps the most dramatic feature in Fig. [2] is the for- 
mation of a bifurcation in the congruence of geodesies. 
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As is evident from Fig. [2] at late times some geodesies 
escape up to the boundary and some plunge deep into 
the bulk. Separating escaping from plunging geodesies is 
precisely one geodesic that does neither. This geodesic, 
shown as the solid blue curve in the figure, defines the 
location of a null surface beyond which all events arc 
causally disconnected from observers on the boundary. 
This surface is the event horizon of the geometry. 

After the time r/ = 2.25, the boundary geometry be- 
comes flat and unchanging, no additional gravitational 
radiation is produced, and the bulk geometry approaches 
a slowly evolving form. The rapid relaxation of high fre- 
quency modes can clearly be seen in the behavior of A/r 2 
shown in Fig. [2] — all of the high frequency structure 
in the plot appears only during the time interval where 
the boundary geometry is changing and creating grav- 
itational radiation. Physically, the rapid relaxation of 
high frequency modes occurs because the horizon acts as 
an absorber of gravitational radiation and low frequency 
modes simply take more time to fall into the horizon than 
high frequency modes. Therefore, as time progresses the 
geometry relaxes onto a smooth universal form whose 
temporal variations become slower and slower as r — > oo. 

One can systematically construct late-time asymp- 
totic expansions of boost-invariant solutions to Einstein's 
equations [7J. The expansion, which is a power series ex- 
pansion in gradients, is dual to the hydrodynamic expan- 
sion in the field theory. This is natural, as the late time 
evolution of the field theory state in conformal Af = 4 
SYM must be described by hydrodynamics. In the grav- 
itational setting, the metric is expanded in terms of 4c? 
spacetime gradients of slowly varying fields. For the case 
of boost invariant flow, each spacetime derivative intro- 
duces a factor of 1/(At) 2 / 3 into the solution, where A 
is an energy scale which characterizes the initial energy 
density [24]. The numerical coefficients of the expansion 
are related to transport coefficients in the dual gauge the- 
ory, and are independent of the initial conditions used to 
create the black hole geometry. Therefore, at asymptoti- 
cally late times all sensitivity to the details of the initial 
conditions used to created the black hole geometry is 
isolated within the energy scale A, up to exponentially 
decreasing corrections to the late time behavior. 

At asymptotically late times, the boost invariant gra- 
dient expansion of Ref. 24J yields a metric 



ds 2 



-11— 



dT 2 



dx 



2drdr, (21) 



where tvi(t) « ttA/(At) 1 ' 3 is the approximate location of 
the event and apparent horizons, whose positions asymp- 
totically coincide at late times. The asymptotic metric 
(21 1 has a Hawking temperature 



Hawkint 



= A/(Ar) 1 /3, 



(22) 



which is proportional to the horizon radius rft(r). As 
time progresses, the horizon slowly falls deeper into the 
bulk, and the temperature of the black hole decreases 




■Event horizon 
Apparent horizon 



2 4 

T 

FIG. 3: Area of the event horizon and apparent horizon, per 
unit rapidity, as a function of proper time r. The growth of 
the apparent horizon area, shown by the magenta dotted line, 
is causally connected to the changing boundary geometry. In 
contrast, the growth of the event horizon area, shown as a 
solid blue line, is non-zero at arbitrarily early times, long 
before the boundary geometry has started to change. 



as t -1 / 3 . The falling of the horizon into the bulk, as an 
inverse power of r, is clearly visible in the numerical data 
presented in Fig. [2] 

Fig. [3] shows a plot of the area (per unit rapidity) of 
the event and apparent horizons, again for c = 1, as a 
function of t. The area (per unit rapidity) of the appar- 
ent horizon is given by £(r?i(T), t) 3 where r^(r) is the 
apparent horizon location (given by a zero of S). The 
area (per unit rapidity) of the event horizon is also given 
by S 3 , but instead evaluated on the null geodesic defin- 
ing the event horizon. The area of the apparent horizon 
starts off at zero, and grows rapidly for r in the interval 
( T i? r /)- This is to be expected, as it is during this in- 
terval of time that the rapid variation of the boundary 
geometry produces infalling gravitational radiation which 
is subsequently absorbed by the horizon. As radiation is 
absorbed, the horizon area must grow. After the produc- 
tion of radiation ceases, the the geometry relaxes onto 
the asymptotic form (211 and the area (per unit rapid- 
ity) of the apparent and event horizons slowly approach 
a constant. From the figure, one sees that the growth of 
the apparent horizon area changes rather abruptly near 
time Tf. This reflects of the fact that the boundary ge- 
ometry ceases to produce infalling radiation after time 
Tf. The flux of radiation through the horizon decreases 
dramatically after r/ and correspondingly, so does the 
growth of the apparent horizon area. 

In contrast to the apparent horizon area, which is non- 
zero only in the causal future of the boundary time r», the 
event horizon area is non-zero arbitrarily far in the past, 
long before the boundary geometry starts to change. 
This reflects the teleological nature of event horizons. 
The event horizon separates events which are causally 
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■Event horizon 
Apparent horizon 

■2 nd order gradient expansion 



time asymptotic expansions, computed through second 
order in gradients. These asymptotic results are [25l [3TJ 



r 



FIG. 4: Close-up view of the event horizon and apparent 
horizon areas, per unit rapidity, as a function of proper time 
r, together with their corresponding asymptotic expressions 
( |23[ | . Both horizon areas are very well approximated by their 
asymptotic expansions, at second order in gradients, after 
time Tf — 2.25 when the boundary geometry becomes flat. 
Note the rather abrupt change in the growth of the apparent 
horizon area at Tf. 



disconnected from boundary observers. As Fig. [2] clearly 
shows, even before the boundary geometry has started to 
change there are events which are causally disconnected 
from the boundary. These events are, by definition, be- 
hind the event horizon. Simply put, the black hole exists 
before the boundary deformation has begun! 

Because the radial geodesic defining the event hori- 
zon is moving outwards at the speed of light, before the 
boundary geometry starts to change the area of the event 
horizon grows like 4(fc+r) /(fc— r) 3 , where k is the value of 
t + 2/r on the geodesic defining the event horizon. The 
appropriate value of the constant k can only be deter- 
mined when the entire future of the geometry is known. 4 
Because of its acausal nature, the area of the event hori- 
zon cannot correspond to the entropy of the system in a 
non-equilibrium setting. In contrast, it does appear sen- 
sible to regard the apparent horizon area as a measure of 
thermodynamic entropy in a non-equilibrium setting. 

To facilitate a quantitative comparison between our 
numerical solutions to Einstein's equations and the late 
time gradient expansion of Ref. [24], Fig. [4] shows a close- 
up view of the areas (per unit rapidity) of the event and 
apparent horizons, together with the corresponding late- 



4 This manifests itself as follows. At asymptotically late times, the 
location of the event horizon coincides with the zero of A(r,r), 
so the unique outgoing radial geodesic that approaches the zero 
of A(r, 7") as t — > oo defines the event horizon. To locate the 
position of this geodesic at early times, and hence determine the 
horizon area, one must know the entire future of the geometry. 



^EH = ^'A 2 

A AH = ir 3 A 2 



1 - 
1 - 



1 6 + 7r + 61n2 

2tt(At) 2 / 3 + 24tt 2 (At) 4 /3 
1 2 + 7T + In 2 " 



2tt(At) 2 / 3 24tt 2 (At) 4 / 3 



, (23a) 
(23b) 



for the event and apparent horizon areas, respectively, 
up to O ((At) -2 ) corrections. From the figure one sees 
that the asymptotic expansions, shown in the figure as 
the dashed black lines, agree very well with the complete 
numerical results. In fact, at time Tt when the boundary 
geometry becomes flat, the asymptotic forms agree with 
the full numerical results for both horizon areas to within 
0.11%. 

For c = 1, our numerically measured value of A is 
0.8. Consequently the first order corrections appearing 
in Eqs. (23 1 generate 10% corrections at time Tf, while 
the secondorder terms yield 0.20% and 0.56% corrections 
to the event and apparent horizon areas, respectively. 

This comparison shows that the geometry in the bulk 
(as probed by the horizon areas) is already very well ap- 
proximated by the gradient expansion of Ref. |24j at time 
Tf. However, it must be stressed that this very early 
agreement with hydrodynamics is specific to the horizon 
areas, and is not so true of other observables which are 
sensitive to the anisotropy in the geometry, such as the 
SYM stress tensor, which we discuss next. 

We now turn to a discussion of our results for bound- 
ary field theory observables. Fig. [5] shows plots of the 
energy density and transverse and longitudinal pressures 
produced by the changing boundary geometry ([!]), when 
c = ±1. These quantities begin at zero before time Ti, 
when the system is in the vacuum state, and deviate from 
zero once the Ad geometry starts to vary. During the in- 
terval of time where the Ad geometry is changing, the 
energy density generally grows and the pressures rapidly 
oscillate: work is being done on the field theory state. 
After time Tf the boundary geometry becomes flat and 
no longer does any work on the system. As time pro- 
gresses, non-hydrodynamic degrees of freedom relax and 
at late times the evolution of the system is governed by 
hydrodynamics. The late time hydrodynamic behavior 
manifests itself as the smooth tails appearing in Fig. [5j 

The two sets of plots in Fig. [5] contrasting c = +1 and 
— 1, are qualitatively similar, with the main difference 
being the phase of the oscillations in the pressures. For 
example, for c = — 1 the transverse pressure is negative 
at Tf whereas for c = +1 the transverse pressure is pos- 
itive and larger than the longitudinal pressure, which is 
nearly zero at Tf. As local equilibrium requires that the 
transverse and longitudinal pressure be nearly equal |3J, 
one sees that in either case the system is far from equi- 
librium at Tf. Furthermore, from the figure one sees that 
for either sign of c, the transverse pressure approaches 
the longitudinal pressure from above. As we next dis- 
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c=-l 





FIG. 5: Energy density, longitudinal and transverse pressure, all divided by as a function of time for 



-1 (left) 



and c = +1 (right). The energy density and pressures start off at zero at time n = 1/4 when the system is in the vacuum state. 
During the interval of time r £ iji, r f) = (0.25,2.25), the 4d geometry is changing and doing work on the field theory state. 
After time 77 the deformation in the geometry turns off and the field theory state subsequently relaxes onto a hydrodynamic 
description. The smooth tails in both plots occur during this regime. At late times, from top to bottom, the three curves (in 
both plots) correspond to the energy density £, transverse pressure Vx, and longitudinal pressure V\\. 



cuss, this is always the case in the hydrodynamic limit of 
boost invariant flow. 

From the gravitational asymptotic expansion of 
Ref. |24j . one can compute the SYM stress tensor via 
Eq. ((nj). The results read [24] 



3tt 4 A 4 
4(At) 4 / 3 

^ 4 A 4 
4(At) 4 / 3 

^ 4 A 4 
4(At) 4 /3 



1 



1 



2Ci 

(At) 2 / 3 " 

C 2 
3(Ar) 4 /3 

2Ci 



Co 



(At) 4 / 3 



5C 2 



(At) 2 / 3 3(At) 4 / 3 



(24a) 
(24b) 
(24c) 



up to C((At)~ 2 ) corrections. The constant C\ is related 
to the viscosity to entropy density ratio of the plasma, 
while the constant C2 is related to second-order hydrody- 
namic relaxation times. For strongly coupled SYM [25] . 



1 

3^ 



C 2 



2 + In 2 
18tt 2 ' 



(25) 



The form ( 24 ) for the stress-energy can also be obtained 



from hydrodynamic considerations alone, together with 
knowledge of first and second order transport coefficients, 
and the assumption of boost invariance [3H |3S] . 



It is evident from the leading terms of the result ( 24 ) 



that at late times the stress-energy tensor approaches the 
ideal hydrodynamic form 



diag(-3, 1,1,1), 



with a time-dependent temperature 
T(t)=A/(At) x / 3 , 



(26) 



(27) 



which matches the Hawking temperature (22 1 of the 
black brane in the gravitational description. The ideal 
stress tensor (26 1 is completely isotropic. Sublcading 



terms in the result ( 24 1 show that the transverse pressure 
differs from the longitudinal pressure when viscous ef- 
fects are taken into account. In particular, as mentioned 
above, first order viscous corrections make the transverse 
pressure larger than the longitudinal pressure. 

To facilitate a quantitative comparison between our 
numerical results for the stress tensor and the late-time 
hydrodynamic expansions, Fig. [6] shows the energy den- 
sity and pressures for c = 1/4, 1 and 3/2, with the cor- 
responding hydrodynamic forms p4| plotted on top of 
the numerical data. The plots start at time t — Tf. In 
all three plots, one clearly sees the stress-energy com- 
ponents approach their hydrodynamic approximations. 
Moreover, in all plots one sees a substantial anisotropy 
even at late times where a hydrodynamic treatment is 
applicable. In other words, the effect of viscosity is very 
evident in these results. 

Looking at Fig. [6] for time r — Tf and c — 1/4, one 
sees that the transverse and longitudinal pressures are 
almost equal and opposite in magnitude at this time. So 
the system is initially very far-from-equilibrium. How- 
ever, for c = 3/2 the pressures are both positive, and 
system is much closer to equilibrium at Tf. At first sight 
this might seem peculiar: how can it be that for larger 
values of c, where the size of the perturbation in the Ad 
geometry is huge, the system takes less time to reach local 
equilibrium! Qualitatively, this apparent puzzle is easy 
to understand. For large c, the changing geometry does 
more work on the system and consequently the system 
reaches a higher effective temperature. Because SYM is a 



c= 1/4 
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C= 1 
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FIG. 6: Energy density, longitudinal and transverse pressure, all divided by N^/2tv 2 , as a function of time for c = 1/4 (left), 
c = 1 (middle) and c = 3/2 (right). From top to bottom, the curves are energy density (blue), transverse pressure (green), and 
longitudinal pressure (red). The dashed black lines in each plot show the second order viscous hydrodynamic approximation 
( |24[ | to the different stress tensor components. Note the significantly different ordinate ranges in the three plots; the size of the 
difference between the transverse and longitudinal pressure grows with increasing c. 
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1.2 


0.87 


0.89 


1.3 


1.9 
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3.3 
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2.0 


1.7 


1.4 


1.1 


0.84 


0.85 


1.1 


1.5 


1.8 


2.1 


(t*-t/)T* 


0.00 


0.05 


0.11 


0.19 


0.24 


0.24 


0.20 


0.11 


0.04 


0.00 


■P ± (T f )-P U (Tf) 


0.06 


-0.03 


-0.22 


-0.56 


-1.1 


1.6 


0.91 


0.47 


0.24 


0.13 





TABLE I: Quantities characterizing the relaxation to equilibrium, for various values of the boundary perturbation amplitude 
c. The relaxation time r» (in units of A) is the time at which the transverse and longitudinal pressures deviate from their 
hydrodynamic values ( |24[ ) by less than 10%. T* is the temperature at time t*, and A is the scale appearing in the hydrodynamic 
expansion (24 1 (both measured in units of A" 1 ). The quantity (r»— Ti)T* measures the total time in units of T, required to 



produce the plasma and relax to near local-equilibrium. The quantity (t* — T/) Tl* measures the time in units of T, required 
for the plasma to relax after the deformation in the geometry ceases. The quantity [P±(Tf) — P|i(t/)1 /£(jf) is the pressure 
anisotropy, relative to the energy density, at time Tf. 



conformal theory, relaxation times for non- hydrodynamic 
degrees of freedom must scale inversely with the temper- 
ature, and hence must vanish as the local energy density 
diverges. Therefore, in the c — > oo limit the system will 
always be very close to local equilibrium — even while 
the Ad geometry is changing — and the anisotropy in 
the pressures will vanish immediately at Tf. As a con- 
sequence, one learns little about the physics of the re- 
laxation of non-hydrodynamic degrees of freedom in the 
c — > oo limit. 

Table |l]shows how various quantities characterizing the 
relaxation of the plasma depend on the boundary pertur- 
bation amplitude c, within the range [—2,2]. Included in 
the table is the time t*, beyond which the stress ten- 
sor agrees with the hydrodynamic approximation ( 24 1 to 
within 10%. Also shown is the temperatures J 1 * at time 
t*. the scale A measured in units of t„, and the time 



intervals r*— Tj and T*—Tf measured in units of 

From the table, one sees that as the magnitude of c 
increases, so does the temperature T*. Moreover, as the 
magnitude of c increases, one sees that the time scale r* 
approaches t/ = 2.25. In particular, for |c| = 2 the stress 
tensor is already within 10% of its hydrodynamic limit 
at Tf. As discussed above, both of these features are 
to be expected. Increasing |c| means that the changing 
geometry does more work on the system, producing a 
larger energy density, and consequently the relaxation 
times of non-hydrodynamic degrees of freedom decrease. 
In all cases presented in Table [I] the relevant dynamics 
— from the production of the plasma to its relaxation to 
near local equilibrium (where hydrodynamics applies) — 
occur over a time t* — r, < 2/T*. 

From Table [!] one also sees that for |c| < 1/2 the time 
scale r* at which a hydrodynamic treatment becomes ac- 
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curate occurs when At* sa 1. For larger values of |c|, 
At* is bigger. Examining the size of the coefficients in 
Eq. p4| shows that the hydrodynamic expansion is quite 
well-behaved for At > 1. From the time-temperature re- 
lation (27 1, one can convert At* > 1 to the estimate 
t*T* > 1 which, from Table [fl, is indeed the case. 



Examining the size of the cofficicnts in the scries ( 24 ) 
shows that the second-order (At) -4 / 3 terms are quite 
small compared to the leading (At) -2 / 3 viscous terms 
when At > 1; they only become comparable when At 
is below 0.1. Hence, the fact that hydrodynamics is 
not accurate until At is larger than 1-2 (depending on 
the value of c) indicates that the physics which deter- 
mines the onset of hydrodynamic behavior is not respon- 
sible for higher order terms in the hydrodynamic expan- 
sion becoming comparable to lower order terms. Rather, 
the change in behavior from non-hydrodynamic far-from- 
equilibrium behavior to near-local-equilibrium hydrody- 
namic response must be reflecting the relative importance 
of exponentially relaxing non-hydrodynamic degrees of 
freedom in comparison to the slowly relaxing hydrody- 
namic modes. This means that one cannot accurately 
identify the domain of utility of the hydrodynamic de- 
scription by asking when the late-time gradient expan- 
sion breaks down. A similar conclusion was also reached 
in Ref. |36j by analyzing small perturbations on top of 
an infinite static plasma. 

It is instructive to discuss the qualitative origin of the 
relaxation time t* from the perspective of the 5d gravi- 
tational problem. First consider the limit |c| — > oo. In 
this limit, large amounts of gravitational radiation are 
produced by the changing boundary geometry and the 
amount of energy which falls deep into the bulk diverges. 
As a consequence, the horizon radius must approach the 
boundary as |c| — » oo. The infall time for radiation to 
travel from the boundary to the horizon is roughly equal 
to the inverse horizon radius (in our coordinate system) , 
so in the \c\ — > oo limit, the geometry outside the hori- 
zon effectively responds instantaneously to the changing 
boundary geometry. Therefore, as |c| — > oo, the system 
needs no to time to return to local equilibrium after the 
geometry stops changing at time t/. 

Now consider the \c\ — > limit. For small |c|, dis- 
tinct dynamics occurs on the time scales A and A / yjcf 
[with A = \{tj — Ti)]. First of all, irrespective of how 
small c is, the positions of the apparent and event hori- 
zons are rapidly varying only over the time scale A. This 
is because it is during the time interval Ti < t < tj 
that gravitational radiation is being produced and ab- 
sorbed by the horizon, creating most growth in horizon 
area. The parametric size of the horizon radii at t/ is 
~ y/\c\/A. Qualitatively, this makes sense since little 
radiation is produced and very little radiation falls into 
the bulk of the geometry when |c| is small. Hence, the 
black hole size will vanish as |c| — * 0. An easy way to 
understand the y/\c\ scaling is to note that the total en- 
ergy added to the field theory state cannot depend on 
the sign of c, and therefore must be quadratic in c in the 



small |c| limit. The final state energy density in the field 
theory will scale as r\, so the horizon radii must be pro- 
portional to y+c|. Consequently, after time Tf it takes 
a time ~ A/-^/[c[ for any remaining short wavelength 
perturbations to fall into the horizon. It is during this 
interval of time that the geometry undergoes its relax- 
ation onto the slowly evolving hydrodynamic form. We 
therefore see that t*T* should have a non-zero 0(1) limit 
when \c\ — ► 0. 5 

Last, we discuss the relevance of our work to more com- 
plicated numerical relativity problems in gauge/gravity 
duality. As discussed in the Introduction, an interest- 
ing future direction is the study of collisions of gravita- 
tional shock waves in AdSs, as this is dual to the collision 
of sheets of matter in SYM and mimics the collision of 
large, highly boosted nuclei in heavy ion collisions. In 
the simplest setting, one can study shock waves which 
are translationally invariant in two transverse directions 
|16j . The corresponding gravitational problem is there- 
fore 2+1 dimensional. While we have studied a sim- 
pler 1 + 1 dimensional gravitational problem in this pa- 
per, there are several lessons which may provide insight 
relevant for more difficult problems. First, we found it 
necessary to solve Einstein's equations analytically near 
the boundary with a power series expansion in the ra- 
dial coordinate. This was required as the presence of the 
negative cosmological constant makes the near-boundary 
geometry singular. More specifically, careful asymptotic 
near-boundary analysis was required to determine the 
appropriate subtraction terms needed to make the inte- 
grals (13) finite and produce a numerical scheme which 



remains accurate near the boundary. The same issue will 
arise in gravitational problems with less symmetry. 

Another important lesson concerns the choice of coor- 
dinates. Because of the presence of the negative cos- 
mological constant, all matter and radiation tends to 
fall inward into the bulk. This universality of gravita- 
tional infall motivates the use of coordinates specifically 
adapted for inf ailing motion. The generalized infalling 
Eddington-Finkelstein coordinates we used, which assign 
a constant "time" coordinate to all events on infalling 
null radial geodesies, are especially appropriate for nu- 
merics. Had we used a time coordinate which defined 
a spacelike slicing of the geometry, then we would have 
wasted computational time solving for the geometry deep 
in the bulk before any signals from the boundary had ar- 



A brief comment on the relation between our work and the re- 
cent paper of Bhattacharyya and Minwalla |22| may be in order. 
These authors examined black hole formation and thermaliza- 
tion in response to an arbitrarily weak boundary perturbation 
coupling to the dilaton. A noteworthy finding in this work was 
"instant thermalization" (as probed by measurements of local 
operators) after the boundary perturbation turned off. However, 
this is the case for asymptotically AdS4 spacetime and, as clearly 
stated in Ref. |22| . is not expected to hold in asymptotically AdSs 
spacetime or for non-infinitesimal boundary perturbations. 
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rived. With a null time coordinate, a signal propagating 
in from the boundary at r = oo arrives "instantaneously" 
at r = 0. Moreover, the generalized inf ailing Eddington- 
Finkclstein coordinates yield a metric ^ which is non- 
singular on the horizon. Coordinates which do not yield a 
metric regular at the horizon, such as Fefferman-Graham 
coordinates, are not well suited to numerical initial value 
problems. 

V. CONCLUSIONS 

Using gauge/gravity duality, we have studied the pro- 
duction and relaxation of a boost invariant plasma in 
strongly coupled Af = 4 supersymmetric Yang-Mills the- 
ory The production mechanism is a time-dependent 
deformation of the four dimensional geometry in which 
the field theory lives. The deformation, which was con- 
fined to a compact interval of proper time, does work on 
the system and thus excites the initial state, which we 
took to the JV = 4 SYM vacuum. Within the context of 
gauge/gravity duality, this problem maps into the prob- 
lem of black hole formation in five dimensions. By solv- 
ing the corresponding gravitational problem numerically, 
and using the gauge/gravity dictionary, we were able to 
compute the field theory stress tensor at all times, from 
the first excitation of the initial vacuum state to the late- 
time onset of hydrodynamics. We found that the entire 



process of plasma creation — from the initial vacuum 
state to the relaxation onto a hydrodynamic description 
— can occur in times as short as one to two times 1/T*, 
where T* is the local temperature at the onset of the 
hydrodynamic regime. We also demonstrated that the 
time at which a hydrodynamic treatment first becomes 
valid does not coincide with the point where the hydrody- 
namic gradient expansion breaks down. This reflects the 
fact that in a far-from-equilibrium state there are non- 
hydrodynamic degrees of freedom. These modes relax 
exponentially, and their relative importance determines 
the onset of the hydrodynamic regime. 

This work, together with our earlier paper |21j . pro- 
vide novel additions to the very sparse set of examples 
of genuinely far-from-equilibrium processes in quantum 
field theory which can be studied with complete theoret- 
ical control. Using techniques similar to those presented 
in this paper, it should be possible to study more de- 
manding problems which have less symmetry. 



Acknowledgments 

This work was supported in part by the U.S. Depart- 
ment of Energy under Grant No. DE-FG02-96ER40956. 
We are grateful to Michal Heller, Andreas Karch, and 
Paul Romatschkc for useful discussions. 



[1] P. Arnold, G. D. Moore, and L. G. Yaffe, "Effective 
kinetic theory for high temperature gauge theories," 
JHEP 01 (2003) 030, |arXiv :hep~ -ph/0209353] 
[2] P. Arnold, G. D. Moore, and L. G. Yaffe, "Transport 
coefficients in high temperature gauge theories. II: 
Beyond leading log," JH EP 05 (2003) 051, 
arXiv : hep-ph/0302165 
[3] P. Arnold, J. Lenaghan, G. D. Moore, and L. G. Yaffe, 
"Apparent thermalization due to plasma instabilities in 
quark gluon plasma," \Phys. Rev. Le tt. 94 (2005) | 
| 072302] |arXiv : nucl-th/0409068| 
[4] P. Arnold, G. D. Moore, and L. G. Yaffe, "Transport 
coefficients in high temperature gauge theories: (I) 
Leading-log results," JH EP 11 (2000) 001, 
|arXiv:hep-ph/0010177| 

S. Jeon and L. G. Yaffe, "From quantum field theory to 
hydrodynamics: Transport coefficients and effective 
kinetic theory," \Phys. Rev. D53 (1996) 5799-5809 
arXiv : hep-ph/9512263 

P. Kovtun, D. T. Son, and A. O. Starinets, "Viscosity 
in strongly interacting quantum field theories from 
black hole physics," \Phys. Rev. Lett. 94 (2005) 111601] 



[5] 



[6] 



[ 



arXiv:hep-th/0405231 

[7] S. Bhattacharyya, V. E. Hubeny, S. Minwalla, and 
M. Rangamani, "Nonlinear fluid dynamics from 
gravity," | JHEP 02 (2008) 045| |arXiv: 0712 .2456 



[hep-th] 

[8] P. M. Chesler and L. G. Yaffe, "The stress-energy 
tensor of a quark moving through a strongly-coupled 



N — 4 supersymmetric Yang-Mills plasma: comparing 
hydrodynamics and AdS/CFT," \Phys. Rev. D78 (2008)| 
0450131 [arXiv : 0712 . 0050 [hep-th] | 
[9] P. K. Kovtun and A. O. Starinets, "Quasinormal modes 
and holography," \Phys Rev. D72 (2005) 0860091 
|arXiv:hep r th /0506184 
[10] J. M. Maldacena, "The large N limit of superconformal 
field theories and supergravity," Adv. Theor. Math. 
Phys. 2 (1998) 231-252, |frrXiv:hep-th/9711200| 
[11] E. Witten, "Anti-de Sitter space and holography," Adv. 
Theor. Math. Phys. 2 (1 998) 253-291, 
|arXiv:hep r th /9802150 
[12] S. S. Gubser, L R. Klebanov, and A. M. Polyakov, 
"Gauge theory correlators from non-critical string 
theory," Phys. Le tt. B428 (1998) 105-114, 
|hep-th79 802109 
[13] E. Shuryak, "Why does the quark gluon plasma at 
RHIC behave as a nearly ideal fluid?," Prog. Part. 
Nucl. Phys. 53 (2004) 273-303, |hep-ph/0312227| 
[14] E. V. Shuryak, "What RHIC experiments and theory 
tell us about properties of quark-gluon plasma?," Nucl. 
Phys. A7 50 (2005 ) 64-83] |arXiv:hep-ph/ 0405066 



[15] U. W. Heinz, "T hermaliza tiQiTat RHIC," \AIP Conf. 

Proc. 739 (2005) 163-180||arfiv:nucl-th/0407067| 
[16] D. Grumiller and P. Romatschke, "On the collision of 
two shock waves in AdS 5 ," \JHEP 08 (2008) 027| 
|arXiv: 0803 .3226 [hep-th] 
[17] S. S. Gubser, S. S. Pufu, and A. Yarom, "Entropy 

production in collisions of gravitational shock waves and 



12 



of heavy ions," \Phys. Rev. D7 8 (2008) 066014 
|arXiv: 0805. 1551 [hep-th]"[ ~ 

[18] L. Alvarez-Gaume, C. Gomez, A. S. Vera, A. Tavanfar, 
and M. A. Vazquez-Mozo, "Critical formation of 
trapped surfaces in the collision of gravitational shock 
waves," |arXiv: 0811 .3969 [hep-th]| 

[19] S. Lin and E. Shuryak, "Grazing collisions of 

gravitational shock waves and entropy production in 
heavy ion collision," |arXiv: 0902. 1508 [hep-th]| 

[20] S. S. Gubser, S. S. Pufu, and A. Yarom, "Off-center 
collisions in AdSs with applications to multiplicity 
estimates in heavy-ion collisions," arXiv : 0902 .4062 

| [hep-th] | 

[21] P. M. Chesler and L. G. Yaffe, "Horizon formation and 
far-from-equilibrium isotropization in supersymmetric 
Yang-Mills plasma," |arXiv: 0812 .2053 [hep-th] | 
[22] S. Bhattacharyya and S. Minwalla, "Weak field black 
hole formation in asymptotically AdS spacetimes," 
|arXiv: 0904 .0464 [hep-th]] 
[23] N. D. Birrell and P. C. W. Davies, "Quantum Fields in 

Curved Space,". Cambridge Univ. Pr., UK. 
[24] R. A. Janik and R. B. Peschanski, "Asymptotic perfect 
fluid dynamics as a consequence of AdS/CFT," Phys. 



strongly c oupled N=4 plasma," | JHEP (2008) 133 
|arXiv: 0805. 2570 [hep-th] | 



Rev. D73 (2006) 045013, arXiv:hep-th/0512162 

S. Kinoshita, S. Mukohyama, S. Nakamura, and K.-y. 
Oda, "A holographic dual of Bjorken flow," Prog. Theor. 



Phys. 121 (2009) 121-164, arXiv : 0807 . 3797 [hep-th] 



[26] M. P. Heller, P. Surowka, R. Lo ganayagam, 
M. Spalinski, and S. E. Vazquez, "Consistent 
holographic description of boost-invariant plasma," 
Phys. Rev. Lett. 102 (2009) 041601||arXiv: 0805 .3774 | 
[hep-th] 

[27] G. Beuf, M. P. Heller. R. A. Janik, and R. Peschanski. 
"Boost-invariant early time dynamics from AdS/CFT," 
|arXiv: 0906. 4423 [hep-th]] 
[28] S. de Haro, S. N. Solodukhin, and K. Skenderis, 
"Holographic reconstruction of spacetime and 
renormalization in the AdS/CFT correspondence," 
Commun. Math. Phys. 217 (2001) 595-622[ 
arXiv : hep-th/000 2230 
[29] P. Anninos, G. Daues, J. Masso, E. Seidel, and W.-M. 
Suen, "Horizon boundary condition for black hole 
space-times," \Phys. Rev. T>51 (1995) 5562-5578[ 
arXi v : gr-qc/ 9412069 
[30] R. M. Wald, "General relativity," . Chicago, Usa: Univ. 
Pr. 491p. 

[31] S. Nakamura and S.-J. Sin, "A holographic dual of 
hydrodynamics," JHEP 09 (2006) 020, 
arXiv : hep-th/0607123 
[32] P. Figueras, V. E. Hubeny, M. Rangamani, and S. F. 
Ross, "Dynamical black holes and expanding plasmas," 
|arXiv: 0902 .4696 [hep-th]] 
[33] I. Booth, M. P. Heller, and M. Spalinski, "Gravity dual 
to boost-invariant flow as a slowly-evolving geometry," . 
To appear shortly. 
[34] J. D. Bjorken, "Highly relativistic nucleus-nucleus 

collisions: The central rapidity region," \Phys. Rev. D27] 



I (1983) 140-151| 

[35] R. Baier, P. Romatschke, D. T. Son, A. O. Starinets, 
and M. A. Stephanov, "Relativistic viscous 
hydrodynamics, conformal invariance, and holography/ 
\JHEP 04 (2008) i00] |arXiv : 0712 . 2451 [hep-th] | 

[36] I. Amado, C. Hoyos-Badajoz, K. Landsteiner, and 
S. Montero, "Hydrodynamics and beyond in the 



